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We refute the criticisms of our work on strong-coupling in the presence of an incoherent pumping. 



Our description of strong-coupling in semiconduc- 
tors [TJ HJ |3] has been recently criticized by two groups: 
Ridolfo et al. jl] and Yao et al. jS] . 

Both groups of authors make the same statements: 
they claim that our master equation is flawed, on the 
ground that its domain of convergence is bounded, they 
both propose to use exclusively a thermal bath for the 
reservoir of excitations of the cavity instead of our most 
general case (Yao et al. require a thermal bath also for 
the excitonic reservoir whereas Ridolfo et al. allow inde- 
pendent pumping and decay coefficients for the Quantum 
Dot (QD)) and they all claim a better agreement than 
our model with experimental data. 

We show in this comment that their criticisms are er- 
roneous and that the alternatives they propose to fulfill 
them are well known particular cases that, in their ap- 
proximations, are also erroneous. 

Both groups focus on our boson model only, and in fact 
only on our Letter on the topic pQ. Many of their state- 
ments are already addressed in our full boson text [2] and 
more crucially in its fermion counterpart [5]. Through- 
out, we made clear that the boson model is adequate 
either in the limit of small pumping, or in the limit of 
bosonic 0D system. The interest of the boson case is in its 
analytical solutions, that allow to explain transparently 
fundamental features of pumping, such as the effect of 
the effective quantum state on the spectral shapes. We 
will show that the analytical solutions for the fermion 
models from [1] and [5] are in fact valid only in trivial 
cases (namely, the boson limit or the uncoupled limit). 

In the following, we shall clarify, on the one hand, that 
our approach does not suffer from any inconsistency or 
pathology, and on the other hand, that in the particu- 
lar cases where the reservoirs of excitations are thermal 
baths, our model also applies (by enforcing the ad-hoc 
constrain in the equation) and show the errors made by 
the approximations of Refs. |H[S]. 



I. VALIDITY OF OUR MASTER EQUATION 

Our model describes the linear regime (of vanishing 
excitations) analytically and the nonlinear regime semi- 



analytically. Its master equation reads: 
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where a is the cavity mode and a the exciton in the QD. 
The cavity is always a bosonic mode. In Ref. pQ and 0, 
a is also a Bose operator, that we note b for clarity, while 
in Ref. [3 , a is a fermion operator, describing a two- level 
system [M] . This choice of a as a Bose operator addresses 
two important cases: i) the limit of vanishing pumpings 
(even if the QD is indeed a fermion emitter) and ii) the 
case where the excitons follow bose statistics. The lat- 
ter case could be realized in large QDs, that recover the 
physics of quantum wells where excitons are known to 
behave as good bosons [BJ. 

The main reproach of [U |S] is that in Eq. (JTJ, the 
effective pumping rates P a ,a can vary independently of 
the effective decay rates 7 a . CT [35 . They observe that this 
can lead to some divergences beyond some critical values 
of pumping, as we have ourselves discussed before, and 
they conclude that the model is flawed. They propose 
instead to use thermal reservoirs, that do not exhibit 
such divergences, at all values of pumping. 

A thermal reservoir for a bosonic mode a (with fre- 
quency cu a and H a = Lo a a* a) at temperature T leads to 
the effective rate of excitation: 



P a = K a tl T . 



(2) 



with ut given by the reservoir Bose-Einstcin distribution. 
It vanishes at T — 0. In thermal equilibrium, the system 
is loosing excitations at a larger rate of: 



7a = K a (l + n T ) = K a + P a . 



(3) 



The parameter n a is the spontaneous emission (SE) rate 
at T = 0. The steady state thermal equilibrium reads 
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At very high temperatures, as the effective income of 
particles approaches the outcome, P a ~ la, the number 
of particles remains finite, since P a < i a . As long as i a ^ 
0, any combination of parameters i a ,P a corresponds to 
a given thermal bath (with n a = "f a — P a and T > 0). 
The lincwidth of the optical spectrum of emission is: 



Ta — la Pa — K a 



(5) 



and is independent of temperature (i.e., of the population 
of the mode), since it is always equal to the spontaneous 
emission decay rate n a . 

It is clear from the above results, that a bosonic ther- 
mal bath cannot provide gain and does not exhibit any 
line-narrowing in its luminescence spectrum. A thermal 
bath is a medium of loss as P a < j a by definition. 

A thermal reservoir is, however, a particular case. In 
out-of-equilibrium conditions, especially under externally 
applied pumping, one can expect deviations from the 
thermal paradigm. 

In contrast to the thermal case, a gain medium can be 
derived with bosonic baths out-of-equilibrium. This is 
discussed in textbooks, e.g., in Chapter 7 of Gardiner and 
Zoller's text [7J. A linear gain can be obtained with an 
"inverted" harmonic oscillator maintained at a negative 
temperature — T' . The effect in the master equation is 
that the effective parameters are now given by a relation 
opposite to Eqs. (2][3): 



j a = G a m^ T ,, P a = G a (l + TO_ T /) = G Q + 7 a . (6) 

In this way, G a is the gain or input of particles into the 
mode at zero temperature. Obviously, given that now 
Pa > "fa, there is no stationary solution to the master 
or rate equations. Since the absence of this stationary 
solution is the source of confusion in Refs |U [S] , we quote 
here Gardiner and Zoller's comment, p°216: 

If 7 < K [that is 7 a < P a with our nota- 
tion], there is no stationary situation, and the 
amplifier gives a signal that increases without 
limit. Essentially, the power being fed into the 
cavity cannot escape fast enough. ( Of course 
the idea of an inverted medium which main- 
tains its inversion independent of power out- 
put is not exactly valid, and depletion effects 
will then need to be considered. The system 
is then essentially a laser). 

Obviously, an ever-growing population will always be 
stopped by some external physical effect (the sample will 
have burnt, the reservoir will be depleted, etc. . . ). There 
is however nothing pathologic in this behaviour. In par- 
ticular, this does not invalidate the results for values be- 
low the critical pumping rates. In the case where i c > P c 
(c = a, b), there is a physical solution to the dynami- 
cal master equation, starting at t = 0: at all (finite) 
times, there is a valid master equation, with positive 
trace, normed to unity, etc. . . However, indeed, the sys- 
tem diverges with time. There is no unphysical behaviour 



or flaw of some sort. Not all dynamical systems have a 
steady state, some because they are oscillatory, others 
because they increase without bounds. 

In our work [TJ ® |3J, we have naturally considered 
the configurations which admit a steady state. We have 
even given analytical solutions for their domain of con- 
vergence, supported by a clear physical picture |36j . 

In a microcavity QED system (a QD in a microcavity) , 
which is a complicated solid state system, open to many 
sources of excitations [El [9], one should also include gain 
effects in the most general case. Granted all together, 
the microscopic coefficients arc likely terms of the form: 



la = «o(l + n r) + G a m_ T , , 
P a = n a n T + G a (l + m_ T /) , 



(7a) 
(7b) 



that is, including loss media and gain media. Net losses 
in the case of a cavity mode comes, among other rea- 
sons, from the fact that the photons can escape the 
cavity through the imperfect mirrors. Net gain could 
come from surrounding off-resonance or weakly coupled 
QDs, high energy QD levels or the wetting layer. A gain 
medium can be obtained by the very configuration real- 
ized with self-assembled QDs in a microcavity. Gardiner 
and Zoller, p°140 (ibid) study a bath of uncorrelated two- 
level emitters that are kept in average in the excited state. 
Quoting them again: 

Note that there is no restriction on and 
N~ [that is, in their notations, the number of 
emitters in the excited and ground states, re- 
spectively, that provide the pump and decay 
terms for the bosonic mode] — this means that 
iV+ > N~ is permissible. Physically such a 
population inversion could only be maintained 
by some kind of pumping, as indeed happens 
in a laser. 

and as indeed could happen in an incoherently pumped 
microcavity. The previous discussion on gain-media that 
applies for the cavity mode can also be extended for 
bosonic emitters (QDs, in our case). In this case, for ex- 
ample, the electron-hole pairs that form excitons inside 
the QD decay from the wetting layer (at higher energy) 
by, e.g., emitting phonons with the corresponding energy 
difference. Such phonons will not be reabsorbed to bring 
back the electron-hole pair to the wetting layer, leading 
to a net source of particles. Net losses take place through 
the spontaneous decay into leaky modes. 

In general, there is thus no reason to restrict the mas- 
ter equation designed to account effectively for the 
largest possible amount of physical effects, to a given 
type of reservoirs (namely, thermal baths). Our study 
aims at the greatest level of applicability and generality 
and provides the tools for the understanding of any case. 
Therefore, better than trying to apply at the theoreti- 
cal level any criteria (other than convergence) to choose 
la, Pa, we prefer to consider them independent. So far, 
the dynamics of lines ( Tb"|lc| ) found its most important 
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domain of applicability with atom lasers and polariton 
lasers [10l [TTJ H21 [13l HH dl [HI [17], that is, systems 
where a condensate (or coherent state) is formed by scat- 
tering of bosons into the final state from another state 
rather than by emission. In both cases, scattering or 
emission, the process is stimulated. In this case the 
income and outcome of particles is a complicated func- 
tion of the distribution of excitons (or polaritons) in the 
higher fc-states. See, e.g., Ref [TT]. In a pulsed experi- 
ment, it is typically time dependent (see, e.g., Fig. 2 of 
Ref. [IS])- Line narrowing is a natural feature of this 
dynamics in such systems. 

It remains, of course, possible that a particular exper- 
iment corresponds to a thermal bath of excitation. In 
this fitting analysis with our model should indi- 

cate this constrain in correlations between the 7 and P 
coefficients. As a result of this analysis, one will then 
understand that the given system refers to the particular 
cases of Eqs. Q and 

The above discussion concerns the bosonic mode. We 
now turn to the fermion mode. 

The thermal equilibrium case, at temperature T, gives 
the counterpart of the boson case (given above): 

la = «<r(l + n T ) = n a + P a , (8a) 
P a = K a n T , (8b) 

where K a is the Einstein ^-coefficient and P a is the Ein- 
stein ^-coefficient. The steady state is the Fermi-Dime 
distribution: 
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The maximum occupation for the emitter is 1/2, at infi- 
nite temperature. It is, therefore, not possible to invert 
the two-level system with a thermal bath, (where, again, 
7 CT > P a ), which is a well known result. In the fermion 
case, even a gain-medium does not lead to a divergence 
thanks to the intrinsic saturation of a two-level system. 
The emission spectrum is also a Lorentzian, with effective 
lincwidth: 



T a = la + P a 



2P„ 



(10) 



which broadens from the decay rate at zero temperature, 
K a when the temperature increases. 

This elementary discussion illustrates the notorious 
fact that thermal reservoirs are unable to achieve pop- 
ulation inversion of a two- level system [19J . Therefore, 
the choices of excitation reservoirs of Yao et al. [5J , for- 
bids lasing in such systems, which is already contradicted 
by experiments |20j . 

The issue of achieving gain with a master equation 
for a fermion emitter has been extensively discussed in 
the lasing literature, in particular in the one-atom laser 
case [HI [221 EU [24J EH EH [271 EH [291 [301 EH [32], and it 
is typically described theoretically also by considering an 
effective negative-temperature thermal bath. A popular 



notation to consider gain and dissipation is [25] : 

P a = T a s 7 7ff = T ff (l - s) . (11) 

With T a > (the linewidth broadening) and < s < 1, 
including both the situations with net losses (s < 1/2) 
and gain (s > 1/2). In this case: 

n* = — P tlT = s. (12) 

la + P<r 

This parameterization is not linked to any particular ex- 
perimental realization but is designed to separate the 
physical effects that lead to line broadening, T a , from 
those that change the population, s. Apart from that, it 
is equivalent to consider directly the effect of varying the 
effective decay and pumping parameters, 7(T and P a , as 
we have done in Refs [3] and [33] . 

The Jaynes-Cummings model couples the fermionic 
and bosonic modes. It has been extensively studied in the 
case of thermal cavity bath and some gain-medium for 
the emitter (which is also the case of Ref. [4]). We studied 
it in its most general form, using the master equation (JTJ . 
Again, if a particular constrain arises from a given real- 
ization of the reservoirs of excitations, such as ([2][3]), pj, 
([7]), Q or some other case, this would appear in our un- 
constrained case from correlations following these trends. 
We expect, as we previously discussed, that a successful 
statistical analysis would indeed inform about underlying 
microscopic details of the excitation scheme. 



II. VALIDITY OF THE PROPOSED 
SUBSTITUTES TO OUR WORK 

In the case of Ridolfo et al. [I] , only the photonic mode 
was excited thermally while the excitonic pump was still 
allowing an unconstrained pumping and decay. Yao et 
al. [5J, on the other hand, require both modes to be ex- 
cited by thermal baths. 

Thermal baths in the linear (boson) model reduce to 
results identical to the spontaneous emission of an ini- 
tial state that is a mixture of excitons and photons in 
the ratio of population n a (t = 0)/nb(t = 0) = P a /Pb- 
The thermal character of the bath merely prevents renor- 
malization of the lincwidths and of the Rabi frequencies. 
The ratio P a /Pb still determines the possibility to resolve 
the line-splitting. This fundamental consequence of the 
effective quantum state is independent of any choice of 
the reservoirs. It is a general result that we have amply 
discussed before [U H] and that is "rediscovered" by 
Ridolfo et al. 

Ridolfo et al. [4 otherwise have mistakes in their for- 
mulas, that certainly bias their analysis. For instance, 
their parameters n a and C should read, in their nota- 
tions: (cf. their Eq. (7) & (8)) 

P a , g 2 (7 

71 = ~t~ 

la la 9 2 (la + lx) 2 + lalx - U> : 
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Also, their spectra are normalized to v2irn a , again ap- 
parently as an error since they compare them directly to 
ours which are normalized to unity. These mistakes do 
not seem to be a misprint, given that the authors state 
in their paper: 

Although at low pump intensities, our ap- 
proach and that of Ref. [Laussy et al.] es- 
sentially represent models of a linear Bose- 
like dynamics of two coupled harmonic oscil- 
lators, nontrivial differences can be appreci- 
ated. 

There should be no difference in the limit of vanishing 
pump. In fact, once corrected as above, their formulas 
and lineshapes do converge to our results at low pumps. 

In the non-vanishing case, of course, the thermal reser- 
voir gives a different result than unconstrained parame- 
ters of our Eq. ([TJ, even in the linear regime (that is, when 
n a -C 1 , although n a is not compulsorily also vanishing) . 
The illustration of this fact was attempted in Fig. (la) 
of Ref. [4], although here also the plot is wrong. With 
non-vanishing cavity pumping, the linear regime can be 
maintained only if the modes are almost uncoupled. This 
is shown in our corrected version (Fig. [IJ. 



Our model (PRB 79:235326, 2009) 
Yao et al. (arXiv:0909.5642) 
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FIG. 1: Cavity emission for the parameters of Fig. 1(a) of 
Ridolfo et al. g|: A = -3.6s, -y a = 1.48s, Pa = 0.49s, 
P a — 0.0078s. Like in their figure — but with the correct 
formulas ( 13 1 — we compare: in solid-blue, the case where 
7a = 1.96s an d i n dashed-purple K a — 1.96s {la = K a + P a , 
thermal bath). Both cases are indeed in the linear regime 
(n a <C 1) since the fermion model [3] gives the same results. 
The system is however almost decoupled due to the large de- 
tuning. The emission is thus basically that of the bare cavity 
(Lorentzian). 

We now turn to the approach of Yao et al. [5J. They 
miss a factor r c in the second term of the denominator of 
their Eq. (4) , although in their case this is a misprint only 
since their figures match the formulas that follow from 
their approximations. These approximations, however, 
are incorrect. 

As we discussed above, it is straightforward in our 
approach to constrain the coefficients, following a given 




FIG. 2: Converged spectra with the full fermion model [3] 
[in solid black] for the choice of pumping reservoirs and pa- 
rameters of Yao et al, along with the approximate spectra 
proposed by these authors [in dashed red] . Their approxima- 
tion is correct only at the smallest values of pumping, where 
it also recovers the boson results of Refs. [U |2] . 



choice of the reservoirs of excitations. If we choose the 
thermal baths advocated by Yao et al, we are able with 
our approach [3] to apply the full fermion model with no 
truncation in the number of excitations, for the parame- 
ters of these authors. We find that their approximation 
is a poor one, as seen in Fig. [2] where we superimpose the 
exact (numerical) result, in solid black, to their approx- 
imate (analytical) formula, in dashed red. As could be 
expected, the agreement is good only at very low pump- 
ing (linear regime) and very high pumping (uncoupled 
regime). It is incorrect in the most relevant region of 
the transition where the doublet collapses (as has been 
reported before, e.g., Fig. 13 of our Ref. [3J). 

Their approximation is also basically flawed in that it 
allows an inversion of population for the two-level sys- 
tem, although it is excited by thermal reservoirs, as seen 
in the magnified version of their Fig. 3(a), that we re- 
produce in dashed lines in our Fig. [3j along with the 
converged solutions from our model [3J (with their pa- 
rameters and choices of reservoirs). Beyond the poor 
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FIG. 3: Converged populations [in solid black] for the cav- 
ity, (a) a) [thin] and the QD, {a + o~~) [thick], for the choice of 
pumping reservoirs and parameters of Yao et al. , along with 
the approximate values proposed by these authors [in dashed 
red]. Again, their approximation is correct only at the small- 
est values of pumping, where it also recovers the boson re- 
sults of Refs. [U |2] . The breakdown of their approximation 
is further manifest by the inversion of population of the QD 
(indicated by the arrow), which is notoriously impossible with 
thermal reservoirs. Our model shows the expected saturation 
bounded by 1/2. 



quantitative agreement when pumping is non- vanishing, 
at the point indicated by the arrow and above, the QD 
is inverted, which indicates a pathology of their approx- 
imation. 

Their implication that cavity pumping is determinant 
to achieve lasing is in contradiction with well-known and 
established facts of the one-atom laser theory. See for in- 
stance our text Ref. |3] where lasing is achieved without 
cavity pumping, thanks to the gain- medium that our gen- 
eral model allows. On the other hand, thermal reservoirs 
of Yao et al. forbids lasing, regardless of the magnitude 
of cavity pumping (note that with their parameters, very 



high cavity populations are already achieved, but they 
have thermal statistics, with second-order correlator </ 2 ) 
that increases rapidly towards 2 with pumping). 

Finally, we want to stress that their Fig. 2(b), that sup- 
posedly represents our model, does not make any mean- 
ingful comparison, since, fitting some data with their 
model [that we have just shown is wrong, but even if 
it was correct], they proceed to plot our boson model 
with their fitting parameters. It is obvious that, the two 
formulas being different, the best-fitting parameters for 
one of them will yield poor agreement on the same data 
for the other. Beside, they should have used our fermion 
model, since they consider a supposedly two-level emit- 
ter in a nonlinear regime. Fitting them independently, on 
the one hand, and comparing them on statistical grounds 
on the other hand, rather than settling for some aesthetic 
of the agreement, is the correct course of action. Fitting 
with the nonlinear fermion model is not a trivial task. 
With E. Cancellieri and A. Gonzalez- Tudela, we have re- 
cently obtained results in this direction, to be published 
shortly. 

In conclusion, we have shown that our work [H [2j [3] 
is correct and that the critics addressed against it [4j |5] 
are unsubstantiated on the one hand, and the proposed 
substitutes are incorrect on the other hand. These au- 
thors do not derive any master equation. They settle for 
thermal reservoirs, which derivation is a standard text- 
book material. This is a particular case of our work that 
they apply incorrectly or beyond its limits of validity. 
We have already provided the valid limit for the nonlin- 
ear regime [3] . There remain many open questions in the 
field. Some can be settled by statistical analysis of exper- 
imental data with our model, Eq. (JT|) , which correlations 
between the fitting parameters can teach about underly- 
ing microscopic mechanisms (such as the nature of the 
bath of excitations, among other). To this intent, we in- 
vite experimentalists to make their raw data available to 
everybody [37] . 
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[36] The limit of validity of the boson master equation is given 
by a total outcome of excitation larger than the income 
7a + 7t > P a + Pb, that is, r a + F b > (in the notations 
of [2]) If one of the oscillators has a net gain (r o < 0), 
there is an extra condition for convergence that has a 
clear physical motive: not only mode a is constrained to 
be neatly lossy, T a > 0, but also the income through the 
"inverted" mode b, given by P„, must be smaller than 
its total outcome, which is given by 7b plus the effective 
(Purcell) decay through cavity emission. That is, Pi, < 

76 + li T ^ where g eS = g/y/1 + [2A/(r„ + F b )] 2 , and 
= g at resonance. 
[37] http://sciencecommons.org/about/towards 



